Phases of Bosons or Fermions in confined optical lattices 



O 

o 

(N 

< 



o 
o 
I 

Oh- 

^ : 



I 

O 
o 



> 

m 
m 
cn 
o 

o 



I 

o 
o 



H. Heiselberg 

(/nil), of S. Denmark, Campusvej 55, DK-5230 Odense M, Denmark 

Phases of Bose or Fermi atoms in optical lattices confined in harmonic traps are studied within the 
Thomas-Fermi approximation. Critical radii and particle number for onset of Mott insulator states 
are calculated and phase diagrams shown in ID, and estimated for 2 and 3D. Methods to observe 
these and novel phases such as d-wave superconductivity are discussed. Specifically the collective 
modes are calculated. 

PACS numbers: 03.75.Ss, 71.10.Pm, 03.75.Lm 



I. INTRODUCTION 

Cold quantum systems in optical lattice have been re- 
alized with Bose and Fermi atoms, and superfluid tran- 
sitions to Mott insulator (MI) solids at filling up to five 
[1-6] . Luttinger and Luther- Emery liquids and MI solids 
appear for ID fermions whereas in 2D it is believed that 
d-wave superconductivity (dSC) appears between metal 
and antiferromagnetic (AFM) phases [7,8], which is re- 
sponsible for the high temperature superconductivity. 
Likewise quantum spin phases [9], supersolids [10], frac- 
tional quantum Hall phases [11,12], mixed [13] and other 
phases may be studied. 

By tuning the barriers of the optical lattices and the in- 
teraction strengths near or far from Feshbach resonances, 
the Fermi [14] and Bose [15] Hubbard models can be re- 
alized in 1, 2 and 3 dimensions [16,17] for any interaction 
and hopping strengths. These are the on-site interaction 
energy U and the hopping parameter t to nearest site. 
Both can be expressed in terms of the lattice potential 
and particle interactions such as the s-wave scattering 
length (see, e.g. [16,17]). The optical lattices are in ex- 
periments usually confined in a trap typically in the form 
of a harmonic potential Vi = V2r^. Here, the lattice point 
ri — id is an integer times the lattice spacing d. The op- 
tical lattices are then clean systems with controlled disor- 
der. A number of MI, superfluid and other phases appear 
as the densities increase towards the center of the trap 
[16-22]. 

It is the purpose here to calculate the phase diagram 
for these models in the Thomas-Fermi approximation 
(TF) at zero temperature, which can be done analytically 
in ID for strong and weak couplings U jt and approxi- 
mated between these limits. The corresponding phases 
in 2D and 3D will also be discussed. Finally, the ID col- 
lective modes are calculated for bosons and fermions in 
harmonically confined optical lattices. 



II. EQUATION OF STATES IN THE 
THOMAS-FERMI LIMIT 



In the limit of many particles in a shallow confining 
harmonic potential, the TF (or local density approxima- 
tion) applies because the length scales over which the 



trapping potential and density varies are long compared 
to phase boundaries [23] and the lattice spacing. As 
will be shown below, the size of the system scales as 
R ~ \/t/V2, and the total number of particles scales 
as ~ R/d. Therefore the average filling fraction scales 
as ^ N/R ~ Ndy/V2/t, and can be held fixed in the TF 
limit: N ^ oo and V2 0. 

In TF the total chemical potential is given by the sum 
of the trap potential and the local chemical potential 
/i(/9) — de/dp, where e and pir) are the local energy and 
density per site, i.e. 



PTot = p{p, U) + V^r" = m(0, U) + 1^2 i?' 



(1) 



It must be constant over the lattice and can therefore 
be set to its value at the edge or radius R of the atomic 
gas, which gives the last equation in (1). As we shall see 
later /x(p = 0, C/) = —zt for both repulsive fermions and 
bosons, where z — 2D is the number of neighbors in D 
dimensions. 



A. ID Fermions 

We start with a ID system of spin-balanced fermions 
with repulsive interactions described by the Hubbard 
model [14] 



H=-tY^ d\^dj. 



u 



(2) 



where, 

At 



is the usual Fermi creation operators, ni„ 



the density, and (ij) sums over nearest neighbors 

The ground state energy per site at zero temperature 
can be calculated [24] from the Bethe ansatz. A good 
approximation for densities p < 1 and all [/ > is [19] 



p{p, U) = -2t cos{tt p//3) . 

I3{U ^ 00) 



(3) 



It is exact in both limits: I3{U 00) — 1 and j3{U — > 
0) = 2. For densities 1 < p < 2 the chemical potential 
is: p{p, U) ^ U - p{2 - p, U). [24] The phase is a band 
insulator (BI) for p = 2, a MI for p = 1 and a Luttinger 
liquid (LL) [25] otherwise. The MI phase has a chemical 
potential gap = U + At cos{n / (3) . It must match [24] 



1 



Aij.{U) = U-4:t + 8t / 

^0 



x[l + cx.p{xU/2t)] 



dx . 



(4) 



where Ji is the 1st Bessel function. Alternatively, one 
could as in [19] require that the energy: e(l, f7) = 
— 2i(/?/7r) sin(7r//3), matches the exact result of the Bethe 
ansatz [24]. We prefer matching the chemical potentials 
because it insures that the Mott transitions are correct. 
The chemical potential is discontinuous at p = 1 by an 
amount 



A/i = [/ - 4t + 8 ln(2)tV'^ + OiU-"^) , 



(5) 



for U ^ t whereas in the weakly repulsive limit the gap 
vanishes exponentially as 



(6) 



A/x = -V[/fexp(-27rt/[/) . 

TT 



The density is now found by inserting (3) in (1). For 
p < 1 we obtain 



p{r) = — arcsin(-\/ R"^ — r'^/Rc 

TT 



(7) 



where Rc = 2^/t/\^. The gap in the chemical potential 
Eq. (1) requires a density plateau at p = 1 [19 22]. We 
can now relate p(0) to the radius R and the corresponding 
total number of atoms 



N 



I p{r)dr 
Jo 



(8) 



in the confined lattice for any U /t. 

We arc particularly interested in the critical radii R^ 
and the corresponding critical number of atoms = 

(2/d) /g " p{r)dr, right when the MI plateaus of filling n 
first appear (-) in the center of the trap and disappear 
(+) forming a MI shell. The MI phase with p(0) = n = 1 
first appear in the center of the trap when N = and 

persist up to iV = where after the core is a LL with 
p > I until TV = N2 where a BI with p = 2 forms in the 
core. 

For U ^ t {f3 ^ 1) the MI first appear in the center 
when Ri — Rc according to Eq. (7), and the correspond- 
ing number of particles is = Nc, where 



4fic _ ^ 
ird ndy V2 



(9) 



Likewise, it follows that for U ^ t, where (3 — 2 and 
A/x = 0, the critical MI radius is shorter R'^ = Rc/V2 
whereas the BI radius is R2 = Rc- The corresponding 
critical number of atoms are 



N+ = (ARc/nd) I aTCCosix^)dx :^0M7Nc, (10) 



and A^2~ = 2iVc 



Returning to U ^ t we must add the chemical poten- 
tial gap A/x = [/■ - 4i to Eq. (1) when p{r) > 1. The 



critical radius, where then density just starts to exceed 
unity p = 1+ in the core is: R^ = Rc\/U/4i, whereas 
the BI with p = 2 appears when i?^ = Rc^/U/U + 1. 
The corresponding number of atoms is 



and 



AT - = Nc 



(11) 



(12) 



We emphasize that the calculated values for the critical 
fillings A^j~ , and and the corresponding radii are 
exact in both the weak and strong coupling limits within 
TF. 

For general U we calculate Ap{U) from Eq. (4) to 
find P{U). Inserting (3) in (1) then gives the density, 
critical radii and the critical fillings as shown in Fig. 1. 
They interpolate smoothly between the analytical results 
U and U ^ f. given above. 

The numerical calculations in Refs. [27,20] find simi- 
lar phase diagrams as the TF continuum limit N ^ 00 
discussed above, however, with a few important differ- 
ences. The finite number of particles makes it impossi- 
ble to resolve the exponentially small gap for ?7 — > 0. 
The MI plateau in the density at p = 1 docs not extend 
over a sufficient number of lattice sites, R^ — Ri ^d, and 
the phase boundary cannot be resolved. Therefore, the 
critical point where the p = 1 MI state vanishes moves 
from [/i = in TF to Ui ~ 0.8t in [20] and Ui ~ 3t 
in [27]. As a consequence the phase diagrams of Refs. 
[27,20] acquire an additional phase transition at J7i and 
N^N^ between phases with and without a "resolvable" 
p = 1 Mott transition. A similar phenomenon will occur 
at finite temperature, where the MI solid melts and the 
critical point moves to a finite U since Ap{U) ~ ksT. 

The attractive Hubbard model (U < 0) tends to 
pair atoms in a Luther-Emery liquid (LE) [26]. The 
trap potential breaks translational invariance forming an 
atomic density wave (ADW) phase with oscillating den- 
sity: pin) = p + Spcos{2Trri/\ADw) [26,19,21]. Here, p 
is the average density, Sp is the amplitude and Xadw = 
2d/p the wavelength. The TF can be applied for the av- 
erage density p as long as R^ )^adw- As described in 
[19] the parametrization 



p. = U/2 - 2a^icos(7rp/2) , 



(13) 



is a good approximation for the EoS at all densities 
and J7 < 0. It satisfies the known properties for weak 
and strong attraction [24], when the prefactor is chosen 

as: a'iU) = tt /~ Jo(.t)Ji (.t)/(1 + cxp(x|t/|/2))((ix/a;). 
It decreases from a{U = 0) = 1 to a{U —00) = 
^7rln(2)t/|f7|. No MI are formed in the LE at p = 1 
whereas a BI appears when p = 2 with radius R2 = 
2at/V2 = aRc- The corresponding number of particles is 



N2 = 2aNc 



(14) 



which is also shown in Fig. 1. 
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FIG. 1. Phase diagram of Fermions in a confined ID optical 
lattice vs. U/t and number of atoms. The phases are denoted 
by the core state. For example, for f7 > the BI phase has 
a band insulator core with p = 2) surrounded by a LL layer 
with 1 < p < 2 again surrounded by a MI layer with p = 1 
and finally a surface LL layer with < p < 1. The curves 
separate the phases where the MI appear (A^i), vanish (A''j^) 
and a BI (iV^) appear in the center (see text). For U <Q the 
LE liquid has a BI core when N > ■ 



B. ID Bosons 

Bosons do not experience Pauli blocking and therefore 
have a rather different phase diagram as fermions oven 
for the same lattice, particle mass and interactions. The 
Bose-Hubbard [15] is given by: 

H = -tY,V^j + \uY,ni{hi-l), (15) 

where bi is the Bose creation operator and Ui = b\bi the 
density. 

Supcrfluid transitions to MI solid phases appear at 
integer filling, n = 1,2,3,..., above a critical coupling 
[27,28] Ui = 3.84t, U2 = 6.2t, and ?7„ = 2.2nt for n > 1. 
Just above a critical coupling a Mott gap opens in the 
chemical potential compatible with a Kostcrlitz-Thouless 
transition: Ay^„ ~ exp{— const /-yU — [/„). When U » 
Un we have A/i„ = U — 2zt, whereas for U < Un there is 
no gap and A/i„ = 0. A good approximation to numeri- 
cal results [27] for U ^ t is the simple parametrization 

H = nU+{U - Ann) sin2[(p - n)7r/2] - zt . (16) 

Here, n = Int[p] is the integer value of the density. The 
EoS is constructed to have the correct chemical potential 
gap A/Lt„ built in as well as the correct critical coupling 
Un- Therefore the EoS and resulting phase diagrams 
are expected to be more accurate than, e.g., calculations 
within the mean field approximation. The compressibil- 
ity dp/d^i vanishes in the MI but diverges for densities 



approaching the MI as in [27] . For U tit has the same 
EoS as that for fermions, namely that of a non- interacting 
gas of a single fcrmion species. The strong on-site repul- 
sion acts as Pauli-blocking similar to the case of strongly 
correlated nuclear liquid at high densities. The EoS can 
give the ID Bose phase diagram at least qualitatively and 
generalizes the one shown in [29] . 

The critical radii for the MI transitions can now be 
calculated in the same way as for fermions by inserting 
Eq. (16) in Eq. (1). In the strong coupling limit {t 0) 
the Boson EoS of Eq. (16) becomes exact and we find 

i?+ = VnU/V2 , (17) 

and 

R- = V[nU-Ann]/V2. (18) 

The corresponding critical number of particles = 

(2/rf) /g " p{r)dr are shown in Fig. 2. When U ^ t 
and 1 

Note, that i?f = Rc and iVf = for U :$> t. 
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FIG. 2. Phase diagram for ID bosons with EoS of Eq. (16). 
The parameter regions numbered n = 0, 1, 2, 3, .. form n MI 
phases of density p = (0), 1. 2, ... n separated by n + 1 metal 
regions between core and surface. In the hatched areas the 
density has a MI core p(0) = 1, 2, 3, .., from below and up. 



C. Higher dimensions 

In higher dimensions D = 2,3 the radii and critical 

fillings may be estimated from the scale of the energy 
and chemical potential. At densities up to unit filling 
the chemical potential varies monotonously from —2zt 
up to or +2zt in the weak and strong coupling lim- 
its respectively. The chemical potential thus varies on a 
scale t and therefore the size of the system remains 



3 



of order _Rj~ ^ ^2ztjV2- Above unit filling the chem- 
ical potential gaps are of order ~ U when U ^ t and 
therefore the scale for the critical radii are unchanged 

R-^^=R+ = ^R^ ,U:^t, (20) 

where i?„ = \/U/V2 is the critical radius for strong cou- 
plings (see Eqs. (17-18)). The corresponding critical 
number of particles is 

1 r " 

-^n^Mj Pi^)d''^ = ^» E ^""^^ ' (21) 

m=0 

where Nu = KniRu/d)^ , with geometrical factor kd = 
2, TT, 47r/3 in D = 1, 2, 3 dimensions respectively. 

The phase diagram may have additional phases as in 
the case of the 2D Hubbard model, which is believed 
to describe high temperature superconductivity. With 
extended interactions and hopping more MI phases at 
fractional filling [12] will appear which may separate into 
supersolid phases [10]. Mixing bosons and fermions leads 
to a number of new phases [13]. Spin imbalance inhibits 
pairing (as isospin asymmetry in atomic nuclei) but may 
also lead to novel pairing phases [32,33]. 

D. High temperature superconductivity 

If high Tc is described by the 2D Hubbard model 
[27,8], then we can expect that the MI is replaced by 
an AFM phase around p — 1 surrounded by two dSC 
phases at densities lower (underdoped) and higher (over- 
doped). According to [30] phases with AFM and dSC or- 
der are mixed when U i^8t but coexist when U^8t with 
a first order transition between two densities. A criti- 
cal point is predicted where a mixed AFM-|-dSC phase 
terminate and is replaced by a first order phase transi- 
tion. Here the AFM and dSC phases are separated and 
coexist with a density discontinuity. In [31], however, 
the pure dSC undergo a first order transition to a mixed 
AFM-|-dSC phase. The densities at which these transi- 
tions take place are aroimd doping x = 1 — p ~ ±0.1. 
If nearest-neighbor interactions or next-nearest-neighbor 
hopping are included, the particle-hole symmetry is bro- 
ken and the phase diagram becomes asymmetric around 
half filling (/o = 1). 

The corresponding phases in confined optical lattices 
will therefore be much more complex than the ID shown 
in Fig. 1. For example, the LL phases will be replaced 
by metal, dSC and AFM phases and mixed phases. Ac- 
cording to [30] there is a also a mixed dSC-l-AFM phase 
between the pure dSC and AFM phases when U^8t 
whereas above the dSC and AFM is separated by a first 
order transition line. The associated density discontinu- 
ity is in stark contrast to the constant density MI phase. 

Whereas the details of the phases around the MI and 
AFM shells are not well known, the surrounding shells 
of dSC phases are again surrounded by two metal phases 



extending to the surface and the BI core respectively. 
It would most interesting to study this phase diagram 
for the 2D Hubbard model on the lattice and to find 
these phases and possible other phases such as spin glass, 
strange metal, striped, checker board, etc., and see how 
they vary with coupling U /t, density and temperature. 

III. OBSERVABLES 

The phases described above can be observed by 
time-of-flight expansion, interference, noise correlations, 
molecule formation, [1-4], collective modes [34,36,37], ro- 
tation [32], etc. We discuss density distributions and 
radii, collective modes and rotation below. 

A. MI shell radii 

The density distribution has recently been measured 

[5,6] for 3-D bosons and the MI plateaus observed directly 
up to ri = 5. The radii of the MI shells were measured 
for different atomic numbers for a large coupling [5]. 

In the strong coupling regime the radii of the MI shells 
n = 1,2, ... are (see Eq. (20) and Ref. [35]) 

Rl=Rl-{n-l)Rl. (22) 

Note that the radii ii„ are continous functions of filling 
whereas the critical radii are related to a critical fill- 
ing iV±. 

The corresponding particle number in D dimensions is 

/ n2 \ D/2 

^ = ^«EU^-(«-l) • (23) 

For bosons the upper limit in the sum is given by n — 1 < 
{Ri/Ruf such that i?„ > 0. For fermions also n < 2. 

The radii and particle number can now be related uni- 
versally, i.e. independent of U , t, V2 and d separately 
as long as we are in the strong coupling regime, when 
scaled by J?„ and Ny, respectively. The curves are shown 
in Fig. 3 for the 3-D case at zero temperature. In 
the recent experiment of [5] a 3D optical lattice is filled 
with atoms in the strong coupling regime such that 
Ru/d = sfufV^^ — 26. The measured radii vs. filling 
is also shown for comparison in Fig. 3 and are in good 
agreement with theory. 



4 




12 3 

(N / 

FIG. 3. Radii R„ for MI shells vs. number of particles for 
a 3D optical lattice. Full curves from left to right show the 
n = 1,2, ..5 shell radii for bosons. The dashed curves show 
the MI n = 1 and BI n = 2 for fermions. The solid, half 
filled and hollow data points are the first three bose MI radii 
measured in Ref. [5]. 



B. Collective modes 

The collective modes has proven very useful observ- 
able in traps [34,36,37] and we will therefore calculate 
the modes for optical lattices. 

The equations of motion are in hydrodynamics and 
for a superfluid given by the equation of continuity and 
the Euler equation. Linearizing around equilibrium these 
equations lead to [23,37] 



+ V(f • v) - f (V-v) , (24) 



for the flow velocity v. Here, the force from the harmonic 

trap potential is f = — V(V2?'^) = —raoj^Y, where we have 
introduced the trap frequency uq through V2 = mujQ/2. 

The collective modes can be calculated analytically 
from Eq. (24) when the EoS is a simple polytrope: 
^{p) — fi{0) (X p'^. At low densities the EoS optical lat- 
tices is on such a form for fermions with 7 = 2/D. In 
ID 7 = 2 for both fermions and bosons according to 
Eq. (3) and (16). At higher densities the eff'ective poly- 
trope 7 = [/'/(/i — p(0))]dp/dp taken at central density 
p{Q) = pc gave rather accurate collective modes in atomic 
traps [36,37] . For ID fermions in an optical lattice we find 
fromEq. (3) 



7 = ^ COt(7r/9e/2/3) : 



(25) 



for Pc < 1. In the weakly interacting and the attractive 
cases U < with EoS given by Eq. (13), we also obtain 
Eq. (25) with P = 2. 

The collective modes of fermions and bosons are the 
same in the strong coupling limit since bosons with EoS 



of Eq. (16) also has the effective polytrope of Eq. (25) 
with (3=1. In contrast a dilute bose EoS {U <^ Ui) has 
7=1 (see Fig. 3). We can expect a continuous change 
of 7 = 1 between these two limits. When pc ^ 1 and 
U ^ t both bosons and fermions have 7 — > whereas 
at Pc = 1 a MI forms in the center making the system 
incompressible (corresponding to 7 ^ 00). 

As in 3D [37] the collective modes can be calculated in 
ID generally for any polytropic EoS. Eq. (24) reduces to 



dx u)q 



l)v = 0, (26) 



where x = r/R. The eigenfunctions arc v cx F{~'n,n + 
27— 1,7-1- l,2x— 1), where F is the hypergeometrical 
function. The corresponding ID eigenvalues for multi- 
poles with n = 0, 1, 2, nodes are 



1 + n[7 + 1 + (n - l)/2] . 



(27) 



The dipole mode (n = 0) has constant flow velocity in 
all dimensions and by insertion in Eq. (24) has frequency 
lo/loq = 1. It is independent of dimension, EoS (7) and 
central density because the dipole mode corresponds to 
the whole cloud sloshing as a rigid body in the trap and 
thus only depends on the trap geometry in one direction 
(as long as pc < 1). The breathing or monopolc mode 
frequency has flow velocity v oc r and by insertion in 
Eq. (24) the monopole frequency is found to be w^/wg = 
2 -|- jD, generally in £)-dimensions. 
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FIG. 4. ID breathing mode frequency uj^ /loq = 2 -|- 7 vs. 
central density. The modes for fermions are labelled with 
U/t = 0, 1, 2, 5, 10, 00. Full curve gives the strong coupling 
limit U ^ Ui for both fermions and bosons. The weak cou- 
pling limit [/ ~ is shown by dashed curve for fermions and 
dotted curve for bosons where ui'^/uii = 3. 

In Fig. 3 we show the ID breathing mode frequency 
(n = 1) ui'^/uiq = 2-1-7 with the eS'ective polytrope 7 
given by Eq. (25). It decreases with increasing density 
until the MI appears in the core which is incompressible. 
Likewise we can expect that the collective modes will re- 
veal the AFM state for the 2D Hubbard model. The dSC 
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may only expected to affect the EoS and collective modes 
weakly. Rotational vortices [32] or RF spectroscopy [38] 
may therefore be better signals for the dSC phases. 



C. Rotation 

By rotating the atoms in the traps the BEC and BCS 
supcrflnids are revealed by vortex formation [32]. Simi- 
larly, for confined 2D optical lattices the layers of BEC 
or dSC superfluids with bosons and fermions respectively 
may be identified by vortices when rotated. The MI 
phase is solid and the atoms in it cannot be rotated unless 
a supersolid phase exists or the lattice itself is rotated. 

The possible density discontinuity [30,31] leads to a 
non-analytic moment of inertia. Due to centrifugal forces 
the critical radii and fillings now also depends on ro- 
tation frequency and therefore the density discontinuity 
can made to occur right in the center by varying either 
the particle number or the angular frequency. Such core 
phase transitions have been considered for neutron stars 
resulting in a characteristic behavior of neutron star spin 
down frequencies with time [39]. 

Solids in confined optical lattices will prefer ellipsoidal 
deformation at fast rotation and varying the angular ve- 
locity may therefore create glitches (MI quakes) as ob- 
served for neutron stars. 



IV. SUMMARY 

In summary, optical lattices provide a clean system 
where 1, 2 and 3D Hubbard and Bose-Hubbard models 
can be realized. The phase transitions between MI, LL, 
LE, superfluids and metals were calculated in ID and in 
2-I-3D for the strong coupling regime, where comparison 
was made to recent data on radii and fillings for bose 
atoms in optical lattices. Collective modes were calcu- 
lated and found to be sensitive to the underlying EoS 
and phases. 2D optical lattices with fermions will there- 
fore give direct insight in the high temperature super- 
conductivity and the dSC-|-AFM phases vs. interaction 
strength, doping and temperature. 
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